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Abstract 

Wave attenuation by ice floes is an important parameter for modelling the Arctic Oceans. At 
present, attenuation coefficients are extracted from linear models as a function of the incident wave 
period and floe thickness. Recent explorations in the Antarctic Mixed Ice Zone (MIZ) revealed a 
further dependence on wave amplitude, suggesting that nonlinear contributions are non-negligible. 
An experimental model for wave attenuation by a single ice floe in a wave flume is here presented. 
Observations are compared with linear predictions based on wave scattering. Results indicate that 
linear models perform well under the effect of gently sloping waves. For more energetic wave fields, 
however, transmitted wave height is normally over predicted. Deviations from linearity appear 
to be related to an enhancement of wave dissipation induced by unaccounted wave-ice interaction 
processes, including the floe over wash. 


1 Introduction 


Waves can penetrate hundreds of kilometres into the ice-covered sea. In doing so, wave energy atten¬ 
uates exponentially with distance of propagation into sea ice. The rate of decay depends on the wave 


period and floe thickness (see, for example, Wadhams et al., 1988 Squire and Moore, 1980). Before 


being completely dissipated, however, waves induce breakup, drift and eventually melt of floes. The 
interaction between waves and ice contributes to the formation and extent of an interface of scattered 
ice at the boundary of open waters and pack ice, generally known as the Marginal Ice Zone (MIZ) 


(Wadhams et al., 1988). The extent of the MIZ plays an important role in modelling Northern and 


Southern Oceans and contributes to global climate. 

In recent years, the seasonal ice retreat has expanded significantly, leaving large parts of the Arctic 


Ocean free (e.g., Perovich et al., 2008 Barber et al., 2009). The emergence of large areas of open water 
in the summer allows longer fetches for wave generation. As a result, increasing significant wave height 
approaches the ice edge (Francis et al., 2011 Asplin et al., |2012' Thomson and Rogers 2014). This 
contributes more substantially to the dynamics of the MIZ (cf. [Kohout et ak 2014) than ever before. 

The majority of waves-in-ice modelling investigations focus on predicting the rate of wave attenua¬ 
tion. The ratio of the prevailing floe diameters to the incident wavelengths determines the mechanisms 
responsible for reducing wave energy. In a field of foes with diameters much smaller than their wave¬ 
length (e.g. pancake ice), wave energy is reduced with penetration distance due to viscous losses (e.g.. 


Wang and Shen, 2011). In contrast, floes with diameters comparable to or larger than the wavelength 


(e.g. floes produced by wave-induced breakup) reflect a proportion of the incident wave energy, dissi¬ 
pate a proportion and transmit the remaining proportion. Under these circumstances, models based 
on the elastic bending of the floes can be used to approximate the underlying physics (see, for exam¬ 


ple, Kohout and Meylan, 2008; Bennetts and Williams, 2010, among others). Models, however, only 
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consider the interaction of the waves with an elastic boundary layer at the surface of the water. For 
simplicity, only small-amplitude waves are considered. Furthermore, other effects such as viscosity, 
floe collisions, submergence of the floe and nonlinear effects are neglected. As a consequence, trans¬ 
mission is simply defined by the proportion of reflected energy, while other sources of dissipation are 
neglected. 

Scattering models have been recently coupled with the spectral wave models to simulate waves in 
Antarctic MIZ (Doble and Bidlot, 2013). Despite limitations, hindcasts are overall consistent with 


buoy data. The model, however, seems to mispredict significant wave height at the peak of the storm 
(see, for example, figure 4 in Doble and Bidlotf 2013). In more severe conditions, the floe is subjected 
to water running over the top of it (overwash), drift, six rigid-body degrees of freedom (i.e. heave, 
surge, sway, pitch roll and yaw), besides elastic motions. These processes enhance wave dissipation 
and their extent is proportional to wave amplitude. Interestingly enough, storm waves recorded in 


the Antarctic MIZ (Meylan et ah, 2014) highlighted how the attenuation rate depends indeed on 


both wave period and wave amplitude (or, in other word, the wave steepness, a proportion of the 
ratio of height to wavelength). It is therefore conjectured that current wave-ice models may misjudge 
attenuation coefficient, especially for large waves. In the present letter, an experimental model of wave 
transmission by an single ice floe is discussed to assess the extent of the departure of the underlying 
physics from model predictions. 


2 Experimental model 

Transmission of an incident wave held by an artihcial ice hoe was monitored in the Extreme Wind- 
Wave Flume of the University of Melbourne, Australia. The hume is 60 m long and 1.8 m wide. Water 
depth was set to H = 0.8 m. The facility is equipped with a cylinder type wave-maker on one end for 
the mechanical generation of waves. A linear beach with slope 1:10 is deployed at the opposite end to 
ensure wave absorption. A detailed analysis of the beach performance revealed a contemned rehected 
waves with energy content up to 10% the energy of the input component. 

A model ice hoe was represented by a h = 10 mm thick polypropylene plastic plate and cut into a 
rectangular shape 1 m long and 1.7 m wide. The Young’s Modulus of polypropylene is 1600 MPa and 
its density is 0.905gcm“^. The plate was deployed at 19.2 m from the wave-maker. No mooring was 
applied to the hoe. The free-hoating conhguration allows drift, overwash, the six rigid-body degrees 
of freedom and elastic motions to occur simultaneously. A schematic of the experimental setup is 
presented in Fig. 

Incident regular waves were obtained by imposing a wave period and amplitude at the wave-maker. 
The wave-maker only produces plain sinusoidal waves. Bound modes typical of Stokes waves appear 
automatically soon after wave generation. Four periods were applied: r = 0.7s, 0.8s, 0.9s and Is. 
These correspond to wavelengths A of 0.76 m, Im, 1.26 m and 1.56 m, respectively. For each period, 
six nominal amplitudes a were imposed. This dehnes six different values of wave steepness ka, where 
k = 27r/A is the wavenumber associated to the wave period. The following values were considered: 
ka = 0.04, 0.06, 0.08, 0.1, 0.12 and 0.14, ranging from gently sloping to storm-like waves. Note that 
wave steepness is a measure of wave nonlinearity and hence it dehnes the intensity of bound modes 
(i.e. second-order-like nonlinearity). 

The water surface elevation was monitored by six capacitance wave gauges at a sampling frequency 
of 1000 Hz. Three of them were deployed in front of the incident edge of the hoe to capture rehected 
wave components. Other three gauges were deployed at a distance of 20.3 m from the incident edge 
(see Fig. to measure the transmitted wave held. This distance ensures enough space for the model 
hoe to drift under the inhuence of waves. From visual observations (based on camera recording and 
markers at the wall), the hoe propagated at a speed between 10 and 50 times greater than the Stokes 
drift, depending on wave amplitude and period. 

The hoe is overwashed by waves. Records of the overwashed water were gathered with two small 
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Figure 1: Schematic of the experimental setup (not to scale). The model floe was represented by a 
polypropylene sheet 1.7 m wide and Im long. No mooring was applied. Waves were generated by a 
cylindrical wave maker at one (right-hand) side and absorbed by a linear (1:10) beach at the other 
(left-hand) side. Water surface elevation was recorded by six capacitance gauges: three in front and 
three in the rear of the floe. Front and rear probes were deployed approximately 20.3 m apart to 
ensure enough space for the model floe to drift under the influence of waves. Photos in the bottom 
panels highlight floe’s behaviour in water. Bottom right panel, particularly, show water running on 
top of the floe (over wash). 
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probes deployed on the top of the floe in a post-experiment phase. For technical reason related to the 
length of wires, the floe was moored during overwash measurements. Overall, overwashed layer was 
observed to range from a few millimetres for gently sloping waves to a few tens of millimetres for the 
storm-like waves. It is worth mentioning that the overwash is generated at the front and rear ends of 
the floes, alternately as they pitched. Steep, shallow-water waves propagated in the overwashed fluid 
itself (see bottom left photo in Fig. [^. When shallow-water waves, travelling in opposite directions, 
up and down the floes, meet, waves often become large enough to break, which dissipates wave energy. 
An example of floe propagation in steep incident waves (r = Is and ka = 0.14) and concurrent 
breaking is shown in the supplementary material. 

For each run, 60s time series of surface elevation were produced. Such a short period of time 
excluded any contamination of reflected waves in the measurement area. Three repetitions of the 
same test were carried out to estimate experimental uncertainties. As the nominal amplitude may 
differ from the one obtained in the flume, the incident wave was recorded at the first 3 gauges in the 
absence of the model floe for benchmarking. 

3 Theoretical model 

Consider a two-dimensional model of a regular incident wave interacting with a floe, with a horizontal 
dimensional and a depth dimension. Let the Cartesian coordinate (x, z) denote locations in the water. 
The coordinate x defines the horizontal location. It points in the direction of the incident wave and 
has its origin set to coincide with the front edge of the floe at rest. The coordinate 2 defines the 
vertical location. It points upwards and has its origin set to coincide with the water surface at rest. 

Linear potential-flow theory is used to model water motions. The water is, therefore, assumed to 
be homogeneous, inviscid, incompressible and in irrotational motion. It follows that the water velocity 
field can be defined as the gradient of a scalar velocity potential, denoted <I>(x, z, t), where t is time. 

Kirchhoff-Love thin-plate theory is used to model the floe, which covers the water surface over the 
interval x G (0, 1). The floe, therefore, bends in response to the wave motion, in addition to responding 
in its rigid-body motions. The rigid motions in the vertical direction consist of translational heave 
and rotational pitch. Further, the floe surges back and forth in the x-direction. Thin-plate theory 
permits the deformation of the floe to be defined in terms of the vertical displacements of its lower 
surface, denoted z = —d + w{x,t), where z = —d is the Archimedean draught of the floe. Surge is 
defined by the horizontal location of the floes centre of mass, denoted u{t). 

Wave amplitudes are assumed to be sufficiently small that linear theory is valid. The water-floe 
system therefore oscillates at the frequency of the incident wave. The velocity potential, displacement 
function and surge can, therefore, be expressed as <h(x, 2 ;,t) = Re{( 5 /ia;)(/>(x, z)e“'‘^*}, w{x,t) = 
Re{C(x)e““^*} and u{t) = Re{^e“'“*}, where iv = I'k/t is angular frequency and cf), C and ^ are 
complex-valued. 

The (reduced) velocity potential, cj), satisfies Laplace’s equation in the water domain, an imper¬ 
meable floor condition and the linearised free-surface condition at points not covered by the floe, 

i.e. 

(pxx + 4>zz = 0 for (x, 2 ) G n, (pz = 0 for 2 ; = —H, 4>z = atf) for x ^ (0, 1) and z = 0, (1) 

respectively, where a = Lo'^/g is a frequency parameter. The water and floe motions are coupled at by 
dynamic and kinematic conditions applied at the wetted surface of the floe at rest. The conditions in 
the horizontal direction are 

fO 

4>x = for x = 0,l and z € {—d,0), and — ahl^ = / dz. (2) 

J-d 

The conditions in the vertical direction are 

(fiz = crC and (1 — (Td)C, -|- FC,'"' = (j) for x G (0, 1) and z = —d, (3) 
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where F = Eh^/{12pg{l — i^^)} is a scaled flexural rigidity of the floe and u = 0.3 is Poisson’s ratio. 

On the incident wave side of the floe, far enough away from the floe that the exponentially decaying 
local motions have died out, the wave held is the sum of the incident, (j)i wave plus a rehected wave, 
4>r- On the opposite side of the hoe, far enough away from the hoe, the wave held is composed of a 
transmitted wave only, (pT- The incident, rehected and transmitted wave potentials are dehned as 


4>i 


cosh kH 




where R and T are the rehected and transmitted amplitudes, which must be obtained as part of the 


solution. Meylan and Squire (1994) devised the hrst method to the solve the above problem, without 
surge, C = 0 and using the shallow-draught approximation, d = 0. Bennetts and Chung (2011) 


developed a solution method for the problem with surge and non-zero draught. The model predicts 
the proportions of incident wave energy rehected and transmitted to be TZ = \R/a\^ and T = |T/ap, 
respectively. The model is conservative, i.e. no energy loss, and thus TZ + T =1. 


4 Results 


Transmitted wave heights were calculated by post-processing records in the lee of the hoe with a 


standard zero-crossing analysis (e.g., Emery and Thomson, 2001). Time series were hrst low- and 


high-pass hltered to remove contaminating components greater than 5.5 times and smaller than 0.35 
times the dominant frequency. Individual, transmitted wave heights were then extracted at each of the 
three gauges and for each of the three realisations. Both down-crossing and up-crossing height were 
considered. To generalise the results, the steepness of individual waves (i.e. ka, where the amplitude 
a is dehned as half the individual wave height H) is considered herein. 

Scatter plots of transmitted versus incident wave steepness are presented in Fig. for the different 
incident wavelengths. A box-and-whisker representation is applied to better express the variability 
of the transmitted hied. The tops and bottoms of each box represents the 25**^ and 75**^ percentiles 
of the samples, respectively. The middle line is the sample median. Whiskers extend to the most 
extreme values. Observations beyond the whiskers (-|- symbols in the hgures) are considered outliers. 
Data uncertainty is found to be clearly affected by the steepness. The largest variability was observed 
for high steep waves, where breaking was most likely to occur. Numerical predictions of transmitted 
steepness, as a function of the incident steepness, are reported too (see dashed line in Fig. [^. We 
remark that model predictions only depends on the wave period/length. 

For the longest wavelength (i.e. X/l = 1.56), the hoe exerts a minimal effect on wave transmission. 


in agreement with held observations (cf. Wadhams et ah, 1988). Overall, experimental observations ht 


model data within the range of experimental uncertainty. The median values (middle red line in the 
box), nevertheless, suggest a weak departure from model predictions, with transmitted wave height 
being slightly overpredicted for steep conditions. Under these circumstances, shallow-water waves on 
the overwashed layer becomes very steep and break, increasing the level of dissipation. Departure from 
model predictions becomes more pronounced with decreasing wavelengths. For X/l = 1.26, the model 
performs well up to an incident ka ~ 0.1. Transmission of steeper waves, however, is substantially 
overpredicted by the models; departure is well outside the most extreme values limits indicated by the 
whiskers. This highlights a substantial dependence on wave amplitude, in agreement with records in 


Antarctic MIZ (Meylan et al., 2014). The model performs least well for X/l < 1.0, where deviations are 
recorded over the entire range of incident steepness. For the shortest waves {X/l = 0.76 m), deviations 
appears to be less significant for incident ka > 0.1. 

We remark that the wave-ice model relies on a balance between reflection and transmission and 
hence other sources of wave dissipation are neglected. To verify that deviations takes place when the 
balance is violated, the sum of a representative reflection and transmission coefficient is presented as a 
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Figure 2: Transmitted versus incident wave steepness: experimental data (box-and-whiskers); model 
predictions (dashed lines). Box-and-whiskers includes all individual transmitted wave steepness for 
the concurrent incident steepness. For each box, the central mark represents the median, while the 
edges of the box are the 25*^ and 75*^ percentiles. The whiskers extend to the most extreme data 
points, which are not considered outliers; outliers are plotted individually (-|- symbols). 
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Figure 3: 


Sum of reflected and transmitted energy versus deviation from linear prediction. 


function of the extent of deviations from model prediction in Fig. Reflection and transmission coef¬ 
ficients are calculated as the change in spectral variance (mO) with respect to the incident conditions 
in front and in the rear of the floe, respectively: 

^ mO^ncI _ .y- y \ttiO rear ^Omc| 

A representative spectrum for each test is obtained by hrst calculating spectra at consecutive windows 
of 4096 points (with no overlap) and then averaging them over the entire time series. Deviations from 
model prediction are estimated as the ratio of median transmitted wave height to the incident wave 
height (/3) for each tests. 

When the extent of transmission coincides with the amount of reflected energy (TZ + T = 1), the 
model performs well and hence no significant deviations are observed (/3 ~ 1). Primarily for steep 
waves, however, this balance is violated (i.e. TZ + T < 1) because of, but necessarily limited to, wave 
breaking on the overwashed fluid. As additional sources of dissipation becomes more relevant, the 
model under performs. 

5 Conclusions 

Experiments in the Extreme Wind-Wave Elume at the University of Melbourne, Australia, were 
conducted to measure wave transmission by a single model ice floe. Observations were compared with 
model predictions to evaluate the extent of model assumptions. The ice floe was modelled with an 
unmoored 10 mm polypropylene plate. Water surface elevation was monitored in front and in the rear 
of the floe. Tests were carried out with regular waves with different periods and amplitudes, ranging 
from gently sloping to storm like waves. 

Analysis of the experimental data indicates that models perform generally well under the effect 
of gently sloping waves. Eor steep (storm-like) wave conditions, additional source of wave dissipation 
like wave breaking on top of overrated fluid affect wave transmission substantially. As breaking 
dissipation violate the assumed balance between reflection and transmission, the model substantially 
under perform. 
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